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A  CLOSED  QUEUEING  NETWORK  MODEL 
FOR  MULTI-ECHELON  REPAIRABLE 
ITEM  PROVISIONING 

by 

Donald  Gross 
Douglas  R.  Miller 
Richard  M.  Soland 


1  1.  INTRODUCTION 

Multi-echelon  inventory  theory  has  been  of  Interest  over  the  last 
two  decades,  both  for  the  theoretical  problems  it  poses  and  for  its 
realism  in  describing  operating  systems.  After  a  relatively  dormant 

period  in  the  late  1960's  and  early  1970's,  a  resurgence  of  interest  j 

occurred.  Most  of  the  recent  work  in  multi-echelon  systems  is  keyed  to  | 

a  model  called  METRIC  that  was  developed  at  the  Rand  Corporation  for  the  j 

U.S.  Air  Force.  This  model  involves  finding  optimal  spares  levels  at  ' 

various  locations  in  a  two-echelon  system,  where  lateral  bases  are  sup-  5 

ported  by  a  single  depot.  More  detail  is  provided  in  the  sections  that 
follow. 

Of  interest  here  is  the  study  of  the  trade-off  possible  among 
spares  levels  and  repair  capacities,  as  well  as  a  more  realistic  model 
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than  is  presently  available  of  the  underlying  stochastic  process  that 
describes  components  which  randomly  fail  and  require  repair. 

The  general  problem  to  be  investigated  is  the  determination  of 
the  optimal  spares  levels  and  repair  capacities  in  a  reparable  item 
multi-echelon  system  in  which  a  finite  number  of  items  is  desired  to  be 
operational  at  any  given  time,  and  in  which  queueing  may  occur  at  the 
repair  facilities  when  all  channels — finite  in  number — are  busy. 

Before  presenting  the  details  of  the  model  developed  here,  we 
first  trace  the  historical  development  of  inventory  theory,  with  partic¬ 
ular  emphasis  on  multi-echelon  efforts,  and  summarize  the  previous  work 
in  repairable  item,  multi-echelon  inventory  control. 

2.  HISTORICAL  PERSPECTIVE 

Inventory  theory  is  said  to  have  begun  with  the  development  by 
Ford  Harris  in  1915  of  the  Economic  Order  Quantity  (EOQ)  model  [see 
HARRIS  (1915)].  The  same  model  was  independently  developed  by  R.  H. 
Wilson  at  about  the  same  time,  and  the  model  is  sometimes  referred  to  as 
the  Wilson  Lot  Size  Formula.  This  simple  deterministic  model  still 
serves  today  as  one  of  the  cornerstones  of  applied  inventory  control. 

In  the  1950's  and  1960's,  interest  in  stochastic  inventory  con¬ 
trol  grew  after  the  publication  of  the  landmark  paper  by  ARROW,  HARRIS, 
and  MARSCHAK  (1951).  A  great  deal  of  more  "sophisticated"  mathematical 
work  then  appeared,  concerned  mainly  with  proving  that  (s,S)  types  of 
control  policies  are  optimal  under  a  wide  range  of  conditions  [see,  for 
example,  ARROW,  KARLIN,  and  SCARF  (1958)].  Most  of  this  work  had  to  do 
with  periodic  review  policies,  that  is,  policies  with  the  decision  rule: 
"When  it  is  time  to  review  inventory,  if  the  inventory  position  (on  hand 
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plus  on  order  minus  backorders)  is  below  s  ,  place  an  order  to  bring 
the  inventory  position  level  up  to  S  Only  a  few  of  the  studies  dur¬ 
ing  that  time  were  concerned  with  how  actually  to  find  optimal  values  of 
the  three  decision  variables  "s,"  "S,"  and  "time  between  reviews." 

In  1959,  GALLIHER,  MORSE,  and  SIMOND  considered  continuous  review 
(s,S)  policies,  whose  decision  rule  is:  "Continuously  monitor  the  inven¬ 
tory  position.  When  it  falls  to  a  level  s  ,  place  an  order  Q  which 
will  bring  the  inventory  position  level  to  S  (Q  ■=  S-s)  ."  These  are 
also  known  as  (r,Q)  models  [see  HADLEY  and  WHIT1N  (1963)]. 

While  there  was  interest  in  multi-echelon  inventory  models  during 
the  late  1950's  and  early  1960's  [see,  for  example,  CLARK  (1958,  1960) 
and  CLARK  and  SCARF  (1960,  1962)],  it  was  not  until  the  1970's  that 
computers  were  able  to  handle  the  difficult  task  of  solving  problems  of 
this  magnitude.  PINKUS  (1971)  extended  the  work  of  Clark  and  Scarf  and 
designed  a  truly  multi-echelon,  multi-product  periodic  review  model  for 
consumable  items,  and  showed  that  "real"  solutions  could  be  obtained. 

The  classic  paper  by  FEENEY  and  SHERBROOKE  (1966)  appeared  during 
this  same  period,  and  ultimately  became  the  basis  of  the  most  popular 
multi-echelon  reparable  item  model  of  today  [see  SHERBROOKE  (1968)]. 

For  reparable  item  control,  a  realistic  model  is  a  special  case  of  the 
continuous  review  (s,S)  policy,  where  s  -  S-l  .  This  is  also  known  as 
a  one-for-one  ordering  policy,  and  is  sometimes  used  in  consumable  item 
inventory  control  for  items  that  are  expensive,  critically  important, 
and  infrequently  demanded.  It  is  a  natural  model  for  reparable  item 
situations  in  that  when  an  item  fails,  it  is  generally  dispatched 
immediately  to  a  repair  facility  and  a  spare,  if  available,  is  "plugged 
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in."  Repairing  the  item  is  analogous  to  ordering  a  new  consumable  item 
from  an  outside  supplier  with  the  repair  time  playing  the  same  role  as 
the  replenishment  leadtime. 

The  METRIC  model  of  Sherbrooke  "multi-echelonized"  the  basic 
(S-1,S)  model  of  Feeney  and  Sherbrooke  by  allowing  a  certain  fraction  of 
the  items  to  be  repaired  at  the  base  and  the  remainder  to  be  sent  to  a 
repair  depot.  The  decision  variables  were  the  levels  of  spares  (the 
S's)  to  be  stocked  in  the  field,  i.e.,  at  each  of  the  bases,  and  at  the 
depot.  MUCKSTADT  (1973)  generalized  METRIC  to  allow  for  a  hieraichical , 
or  indentured,  parts  structure;  the  resulting  model  was  called  MODMETRIC. 

A  key  assumption  of  these  METRIC  models  is  commonly  known  as  the 
ample  service  assumption.  This  means  that  repair  capacity  is  infinite, 
i.e.,  that  there  is  never  any  queueing  of  items  waiting  for  a  repair 
channel.  This  has  the  effect  of  causing  successive  replenishment  lead- 
times  to  be  statistically  independent  and  allows  the  invocation  of  a  pow¬ 
erful  theorem  from  queueing  theory — Palm's  theorem  [see  PALM  (1938)]. 
Palm's  theorem  states  that  as  long  as  there  is  ample  service  (Poisson  or 
compound  Poisson  infinite  calling  population  failure  processes  must  be 
assumed  as  well),  it  is  necessary  to  know  only  the  mean  turn-around  time 
of  failed  items,  and  furthermore,  that  the  steady-state  probability  dis¬ 
tribution  of  the  number  of  units  in  resupply  is  Poisson,  with  parameter 
equal  to  the  mean  number  of  failures  during  an  average  resupply  time — in 
inventory  jargon,  the  mean  demand  over  the  leadtime. 

Other  existing  multi-echelon  repairable  item  models  based'  on  this 
ample  service  assumption  are  ACCLOGTROM  [see  FORRY  (1979)],  SIMPLE  SIMON 
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[see  KRUSE  and  KAPLAN  (1973)],  and  TWOPT  [see  KAPLAN  (1980)].  These 
models  differ  from  each  other  in  their  respective  "bells  and  whistles." 
For  example,  SIMPLE  SIMON  allows  some  old  items  to  be  discarded  and  re¬ 
placed  by  purchases.  ACCL0GTR0M  allows  for  the  modeling  of  reliability 
networks;  that  is,  components  may  be  arranged  in  combinations  of  paral¬ 
lel,  series,  and  k-out-of-n  "circuits."  Some  of  the  models  [METRIC, 
MODMKTRIC,  SESAME  (see  Kaplan,  op.  cit.),  ACCLOGTROM]  also  consider 
finding  optimal  values  of  the  decision  variables,  and  their  mathematical 
optimization  techniques  are  somewhat  dissimilar. 

3.  SUMMARY  OF  PREVIOUS  WORK 

The  previous  models  that  have  the  most  direct  bearing  on  our 
research  are  the  ACCLOGTROM,  METRIC/MODMETRIC,  SESAME,  SIMPLE  SIMON,  and 
TWOPT  models.  SIMPLE  SIMON  and  TWOPT  are  stochastic  process  models 
only;  that  is,  they  give  the  steady  state  probabilities  of  the  numbers 
of  units  in  resupply.  METRIC/MODMETRIC  and  ACCLOGTROM,  in  addition  to 
modeling  the  stochastic  process,  have  methodology  for  finding  the  opti¬ 
mal  spares  levels.  SESAME  basically  concentrates  on  finding  optimal 
spares  levels  and  can  use  METRIC,  TWOPT,  ACCLOGTROM,  or  SIMPLE  SIMON  for 
modeling  the  stochastic  process.  There  are  three  basic  limitations  to 
the  preceding  models: 

•  The  stochastic  process  modeling  is  based  on 

(i)  infinite  source  (calling  population),  and 
(ii)  ample  service  (infinite  number  of  repair  channels) 
assumptions. 

•  Because  of  the  ample  service  assumption,  the  only  decision 
variables  in  these  models  are  spares  levels,  and  thus  no 
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trade-off  considerations  are  possible  among  levels  of  spares 
and  repair  capacities. 

Recent  work  by  GRAVES  and  KEILSON  (1‘17*))  considers  a  single  eeln'lon  sys¬ 
tem  allowing  for  a  general  birth-death  stochastic  process  model,  and  in¬ 
troduces  a  new  set  of  performance  measures  dealing  with  times  for  the 
system  to  recover  after  "failing"  and  times  to  "failure"  when  operating 
satisfactorily.  While  alluding  to  system  design  ramifications,  no  ex¬ 
plicit  optimization  problem  is  formulated. 


4.  PROBLEM  STATEMENT 


The  system  we  study  here  consists  of  a  single  base  (or  group  of 
bases)  with  a  single  base  (or  field)  repair  facility  and  a  single  depot 
repair  facility.  The  problem  can  be  stated  mathematically  as 

Minimize  Z  -  k  y  +  k-c  +  k  c  (1) 

y*cB’cD 

M+y 

subject  to  l  p  >  A  ,  (2) 

n-M  n 

where 


Pn  =  steady-state  probability  that  n  units  are  operational, 

M  =  number  of  components  desired  to  be  operating  (operating 
population  size), 

A  »  minimum  percentage  of  time  all  M  components  are  to  be 
operational  (availability), 
y  =  number  of  spare  components  to  "stock", 

=  base  repair  capacity  in  number  of  channels, 
cl,  ■  depot  repair  capacity  in  number  of  channels. 
kt  *  cost  per  unit  (i«y,D,B)  including  annual  operating  costs  and 

capl  al  investment  amortization  of  a  spare  or  a  repair  channel. 
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The  variables  y,  c„,  and  c~  are  decision  variables  to  be  deter- 
mined  by  an  optimization  algorithm.  The  steady-state  probabilities,  , 

must  be  determined  through  a  stochastic  process  analysis.  We  use  closed 
network  queueing  theory  for  the  latter  and  implicit  enumeration  for  the 
f  ormer . 

5.  STOCHASTIC  PROCESS  MODEL 

The  stochastic  process  can  be  viewed  as  a  network  and  is  shown 
schematically  in  Figure  1. 


Node  U 


Node  B 


Figure  1. — Network  for  a  two-echelon  repairable 
item  system. 


The  network  has  three  nodes,  which  we  refer  to  as  U  ("up"  or 
operating),  B  (base  repair  facility),  and  D  (depot  repair  facility). 
Additional  parameters  a  and  B  are  shown,  where  a  is  the  fraction 
of  failed  items  that  are  diagnosed  as  base  repairable  and  sent  directly 
to  base  repair  (1-ot  is  the  fraction  sent  directly  to  depot  repair).  Of 
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those  that  are  sent  to  base  repair,  a  further  fraction  S  ,  after  under¬ 
going  service,  cannot  be  fixed  and  are  sent  to  depot  repair. 

Holding  tines  at  all  nodes  are  assumed  to  be  independent  exponen¬ 


tially  distributed  random  variables.  At  node  II,  the  holding  time  is  the 
time  to  failure  of  a  component,  with  the  mean  failure  rate  denoted  by 
p  (often  in  queueing  and  reliability  literature,  this  is  denoted  by 
A  ).  At  nodes  B  and  D,  the  holding  times  are  repair  times  and  the  mean 
repair  rates  are  denoted  by  and  ,  respectively. 

Since  all  holding  times  are  exponential,  we  have  a  special  case 
of  a  Jacksonian  network  [see  JACKSON  (1957,  1963)];  this  special  case  is 
a  closed  queueing  network  where  no  items  ever  leave  or  enter  the  system, 
but  circulate  within  the  network  only.  Jackson  (op.  cit .)  showed  for 
general  queueing  networks  with  exponential  holding  times  (which,  because 
of  his  work,  are  often  referred  to  as  Jacksonian  networks),  and  CORDON 
and  NKWELI.  (1967)  showed  for  the  closed  network  case  that  the  joint 
probability  distribution  of  the  number  of  customers  at  each  node  of  the 
network  is  of  product  form.  For  closed  networks,  using  the  notation  of 
BUZEN  (1973),  this  means  that  for  a  k  node  network  with  a  total  of  N 
un i ts , 


»(n1,n2,...,nk)  =  G^N)  ^  jj^)  /  A^n^J  , 


(3) 


where  p(n^,n2, . • . .n^)  is  the  joint  steady  state  probability  that  n^ 
components  are  at  node  1,  at  node  2,  ....  n^  at  node  k  .  The 


x^  are  the  real  positive  solutions  to  the  system  of  equations 


IV 

JjXj  =  PiXiPij  *  j=1’2 . 


(4) 
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*U 


Thus, 


P(nU’nB*riD)  " 


UB 

(l-a+uB)pB 

'  "  V*  * 

XB  1 

*  *D 

auD  • 

1  /PB 

I"  * 

•  1  -f 

no 

G(N)  \ctMy 

/ 

w  L 

.  apD  J  VV  ’ 

where  the  A's  are  given  by  equation  (5),  namely, 

,  n  <  b 


A(n) 


n! 


n!b 


n-b 


n  >  b 


where  b  =  M,  c^,  and  cD  ,  respectively.  Buzen's  algorithm  is  used  to 
calculate  the  constant  G(N);  that  is,  the  value  of  G(N)  so  that 


I  p(ny,nB,nD)  -  1  , 


the  sum  being  taken  over  the  set  S  ,  which  contains  all  triplets 

^nU’nB,niP  SUCh  that  nU  +  nB  +  *D  “  N  *  °nce  t*le  ^oint  Probabilities 

p^nU,nB,nD^  are  obtained*  we  can  calculate  the  marginal  probabilities 

p  ,  which  are  required  for  the  constraint(s) ,  by 
nU 


%  =  \  P^U’VV  * 
s 

where  S1  now  is  the  set  of  all  pairs  (n^.n^)  such  that  n_  +  n^  * 

o  D  B  D 

N-n^  .  This  can  be  done  efficiently  using  Buzen’s  algorithm;  in  fact,  it 
results  as  a  by-product  when  calculating  G(N)  . 

The  above  probability  distribution  is  a  function  of  the  decision 
variables  y,  cD,  and  c.  .  The  distribution  exhibits  certain  monotonic- 
ity  properties  in  relation  to  these  variables;  this  will  play  a  crucial 
role  in  the  optimization  part  of  this  study.  Thus,  before  considering 
optimization,  we  verify  monotonicity. 
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Let  represent  the  number  of  "up"  machines  for  a  steady-state 

network;  the  distribution  of  N  depends  on  (y,c_,c  )  and  thus  we  rep- 

U  D  I) 

resent  the  random  variable  as  N  (y,c_,c_)  .  We  put  the  usual  partial 

li  H  D 

order  on  the  decision  space:  (y,cn,c  )  <  (y*,c’,c’>  if  and  only  if 

y  <  y'  ,  cD<  c'  ,  and  cn  <  c'  .  Now  consider  the  idea  of  stochastic 
JJ  a  L)  - ;  U  % 

ordering  of  random  variables:  N  S<C  N*  if  and  only  if  P(N>  n)  <  P(N'  >  n) 
for  all  n  .  We  can  now  state  the  basic  monotonicity  property  for  the 
steady-state  behavior  of  the  system. 


Theorem.  If  (y.Cg.Cp)  <  (y’.Cg.c^)  ,  then  Nu(y,cB,cD)  <  Nu(y’»Cg,c^)  . 

The  transitivity  of  the  inequalities  implies  that  only  three  cases 
must  be  considered  in  proving  this  theorem:  (i)  (y.c^.cj  <  (y+l.c^.c,.)  ; 
(ii)  (y,cB,cD)  <  (y,cB+l,cD)  ;  and  (iii)  (y,cB,cD>  <  (y,cB>cD+l)  . 

first  consider  case  (i) :  We  must  show  that  the  steady-state  num¬ 
ber  of  "up"  units  Increases  stochastically  when  the  total  number  of  units 
is  increased  by  one  and  the  number  of  repair  channels  is  unchanged.  This 
is  easily  seen  to  be  true  by  modelling  the  system  with  M+y+1  units  as 
a  preemptive  priority  network  with  M+y  high  priority  customers  and  1  low 
priority  customer.  The  distribution  of  the  total  number  of  customers  at 
each  node  will  be  identical  to  the  nonpriority  system  with  M+y+1  cus¬ 
tomers.  The  distribution  of  the  number  of  high  priority  customers  at 
each  node  will  be  identical  to  the  nonpriority  system  with  M+y  customers. 

The  one  low  priority  customer  will  spend  part  of  its  time  at  node  "U"; 

s  t 

thus,  N  (y+l,c0,cn)  a  N  (y,c.,cn)  +  L  ,  where  L  equals  the  number 

U  DU  U  a  U  U  U 

of  low  priority  customers  at  node  "U"  in  steady  state.  Since  >  0 
it  follows  that  Nu(y+l,cB,cD)  ^  Nu*y,CB*cD*  * 
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Now  consider  case  (ii) :  We  must  show  that  the  steady-state  number 
of  "up"  units  increases  stochastically  when  the  number  of  base  repair 
channels  is  increased  by  one  and  all  other  parameters  are  held  fixed. 

This  can  be  demonstrated  by  considering  the  form  of  the  joint  distribu¬ 
tion  of  the  number  of  items  at  each  node,  equation  (3).  Note  that  this 
is  the  conditional  distribution  of  three  Independent  random  variables 
given  that  their  sum  equals  the  total  number  of  items  in  the  system.  Let 
us  denote  these  independent  random  variables  as  Z  (M)  ,  Z_(cD)  and 
Zd(cd)  .  Then,  for  example. 


p<V‘b>  ■  ”> 


*  { 


n 

P0B  n! 


n  i  CB 


(7) 


P0B  - CB  ^ n  ■=  *** 


CB,CB 


where  p0B  is  the  appropriate  normalization  constant.  By  the  product 
form,  equation  (3),  it  follows  that 


St 


Zu(M)  |  Zu(M>  +  ZB(cB)  +  ZD(cD)  =  M+y 


(8) 


It  can  be  shown  by  a  straightforward  algebraic  analysis  using  equation 
(7)  (see  Appendix)  that 


V'b  + 11  ?  w  • 

This  fact  and  equation  (8)  then  imply  that 

Nu(y»cB+1,cD^  »  Nu(y*cB’cD)  ’ 
the  desired  conclusion. 

Case  (iii)  can  be  verified  in  exactly  the  same  way  as  case  (ii), 
completing  the  proof.  The  statement  of  the  theorem  gives  the  desired 
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monotonicity  property  of  the  constraint  function  in  equation  (2).  We  can 
now  proceed  to  the  optimization. 

6.  OPTIMIZATION  PROCEDURE 

The  optimization  aspect  of  this  study  is  a  formidable  one  because 
an  expression  for  the  availability  as  a  function  of  the  decision  vari¬ 
ables  (the  spares  level  and  repair  capacities)  does  not  exist  in  closed 
algebraic  form.  That  is,  the  Pn'8  that  appear  in  the  constraint  (2) 
are  determined  from  the  stochastic  process  model  and  can  only  be  calcu¬ 
lated  numerically  when  the  values  of  y,  c  ,  and  c  are  specified. 

a  U 

The  difficulty  just  indicated,  and  the  fact  that  integer  values 
are  required  for  the  decision  variables,  suggest  the  use  of  an  implicit 
enumeration  scheme  for  the  optimization  algorithm.  One  such  scheme  that 
has  already  been  used  when  closed  algebraic  expressions  were  not  avail¬ 
able  (see  SOLAND  (1973)]  is  that  of  LAWLER  and  BELL  (1966).  However,  it 
requires  that  the  objective  and  constraint  functions  each  be  expressible 
as  the  difference  of  two  monotonic  functions  of  the  decision  variables. 
Thus,  use  of  this  optimization  scheme  interacts  with  the  stochastic  pro¬ 
cess  analysis  in  that  the  latter  is  charged  with  providing  the  required 
monotonicity  properties  of  the  model.  We  have  shown  in  the  stochastic 
process  analysis  that  the  monotonicity  conditions  hold  (the  higher 
y,  cB,  cD  ,  the  greater  the  availability).  The  cost  is  linear  and  there¬ 
fore  monotone.  Thus  we  can  use  the  Lawler- Bell  (L-B)  algorithm.  Use  of 
any  other  optimization  algorithm  would  most  likely  place  similar,  or  even 
more  stringent,  demands  (e.g.,  convexity)  on  the  stochastic  process 
analysis. 
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An  implicit  enumeration  approach  alao  allows  us  to  consider  a 
much  wider  class  of  decision  problems  than  has  heretofore  been  treated 
in  connection  with  multi-echelon  inventory  models.  The  optimization  al¬ 
gorithms  ust.*  in  METRIC,  MODMETRIC,  and  SESAME  are  each  tailored  to  the 
specific  form  of  the  problem  treated,  i.e.,  a  single  specific  constraint, 
either  on  service  level  or  on  budget,  and  are  not  easily  generalized  to 
other  formulations.  Through  use  of  an  Implicit  enumeration  approach, 
however,  we  can  treat  a  variety  of  objective  and  constraint  functions 
and  allow  the  use  of  multiple  constraints.  For  example,  we  can  impose 
the  additional  constraint  of  a  lower  limit  on  the  average  number  of  op¬ 
erating  units: 


M-l 

I 

V1 


+  M 


vM 


>  B  . 


or  a  constraint  on  the  availability  of  a  certain  fraction  of  the  popula¬ 
tion: 


nu' 


Mfy 

l  \ 

,*.9M  "U 


>  A’ 


In  applying  the  algorithm  it  is  necessary  to  have  upper  bounds 
for  the  decision  variables.  Certainly  an  upper  bound  for  both  and 

cD  is  Mfy  (ample  server  case) .  To  get  an  upper  bound  for  y  would 
require  some  knowledge  of  the  particular  application,  for  example,  a 
budget  limit  or  a  manufacturing  or  supply  limit. 

The  algorithm  is  based  on  representing  the  values  of  the  decision 
variables  in  a  single  binary  vector  (a  vector  whose  elements  are  either 
zero  or  one).  Suppose  we  had  a  population  of  ten  items  and  knew  from 
budget  considerations  we  could  afford  at  most  five  spares.  Then  an  upper 
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bound  of  fifteen  would  be  adequate  and  four  bits  (binary  variables)  would 
be  adequate  for  describing  each  variable.  Thus  the  L-B  algorithm  would 
work  with  a  twelve  bit  vector,  which  might  be  arranged  as  follows: 


Hence  in  this  case  the  vector  (0010  |  0010  |  0011)  ,  which  has  value 
2^  +  2^  +  2^  +  tP  =  547  ,  represents  the  solution  y  =  2,  cfi  =  2, 

Cp  =  3  .  The  algorithm  uses  the  binary  vector  whose  value  is  547  in 
determining  which  portions  of  the  solution  space  to  eliminate.  For 
example,  in  the  problem  represented  by  (1)  and  (2),  if  the  preceding 
vector  cannot  satisfy  the  constraints,  no  vector  with  value  less  than  it 
can  either,  and  hence  all  solutions  represented  by  vectors  of  value  less 
than  547  are  eliminated  from  consideration. 

It  is  not  necessary  to  partition  the  vector  into  groups  repre¬ 
senting  each  decision  variable;  y,  cn,  and  c  bits  can  be  intermixed. 

D  U 

For  example,  we  could  use  the  ordering 

^y4*  CB4’  CD4 *  y3’  CB3’  CD3’  y2’  CB2’  CD2’  yl*  CB1*  CD1^ 
where  y^,  c^,  and  cDi  represent  the  ith  bit  of  y,  Cg,  and  cQ  , 
respectively.  If  this  ordering  were  used,  the  vector  with  value  547 


shown  above  would  represent  the  solution  y  *  0010  =  2  ,  c„  =  0001  =  1  , 

O 

Cg  -  1001  *  9  .  Which  ordering  is  most  efficient  to  use  depends  on  the 
problem  and  can  only  be  determined  with  some  experimentation.  The 
reader  is  referred  to  Lawler  and  Bell  (op.  cit.)  for  a  detailed  descrip¬ 
tion  of  the  algorithm  and  further  discussion. 
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7 .  SAMPLE  RESULTS 


The  following  problem  was  run  on  an  IBM  4341  using  Bozen' s  algo¬ 
rithm  to  calculate  the  normalizing  constant  needed  to  yield  the  joint 
probabilities  p(nutnB,nD)  ,  and  using  the  L-B  algorithm  to  find  the 
optimal  solution: 

Minimize  Z  =  20y  +  8cfi  +  10cD 


The  upper  bound  used  on  all  variables  was  2M  and  cases  with  M  “  5, 

10,  20,  and  30  were  solved.  Both  constraints  were  used,  except  for  the 
case  of  M  =  30  ,  where  only  (A^)  was  imposed.  The  results  are  given  in 
Table  1.  Four  different  orderings  were  used  and  the  results  for  the  best 
two  are  shown  In  Table  1,  with  ordering  //I  being  (...,  y^,  c^,  cp^> 

yi-l’  cBi-l ’  coi-i*  and  ordering  if 2  being  (...,  yi>  yi_1 . 

c.,,  cD.  ,,  ...,  c_  . ,  cn,  , ,  ...)  .  Given  in  the  table  are  both  the  CPU 
running  times  in  seconds  and  the  number  of  times  Buzen's  algorithm  was 
required  (number  of  times  the  probabilities  had  to  be  calculated). 

Both  the  L-B  and  Buzen  algorithms  appear  to  be  quite  efficient. 

The  most  demanding  of  the  problems  run  was  the  case  M  =  30  ,  both  be¬ 
cause  it  lias  the  largest  population  size  and  because  it  imposes  only  one 
constraint,  causing  the  L-B  algorithm  to  enumerate  more  solutions  than 
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if  both  constraints  had  been  invoked.  Even  so,  the  problem  took  only 
slightly  over  ten  seconds  to  solve. 


TABLE  1 

SAMPLE  RESULTS 


M 

C*B 

CD 

y* 

z* 

A2 

Ordering  1 

CPU  If 

Ordering  2 
CPU  If 

5 

2 

2 

3 

96 

.938 

.982 

1.45 

41 

0.83 

25 

10 

3 

3 

5 

154 

.926 

.988 

2.97 

64 

1. 35 

38 

20 

4 

& 

•v> 

252 

.907 

.989 

— 

— 

2.55 

66 

30 

6 

8 

11 

348 

.904 

— 

— 

— 

10.11 

137 
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APPENDIX 
PROOF  THAT  :'fi( cg  +  1) 


Even  though  it  is  intuitively  obvious  that  if  one  increases  the 
number  of  servers,  congestion  decreases,  so  that  it  is  certainly  logical 
that  ZB(CB  +  1)  is  stochastically  smaller  than  Zg(cg)  »  it:  is  not 
trivial  to  prove.  We  proceed  as  follows,  using  Equation  (7). 

Consider  the  ratio  P(Zg(cB)  =  n)  /  P{Zg(cg)  =  n  -  l}  ,  which  we 
shall  call  R(cg,n)  .  From  Equation  (7), 


(!b 

I  n 

R(cfi,n)  =  / 

K 

CB 

It  is  clear  that 

R(cfi,n)  >  R(cg  +  1 

This  also  implies  that 

P(Zb(Cb)  -  i)  P<ZB<CB 

p(VcB>  -  Jt  " 


1  £  n  <  cR 

(Al) 

CB  n  5.  M  +  y 


n)  n  =  1,2, . . . ,M+y  .  (A2) 


+  1)  =  1} 

+  1\  ,  TT  0  5  1  <  j  <  M+y  (A3) 


since  ^n_j+j  K(Cg,n)  >  R(cg  +1,  n)  ,  and  hence 

1  ^  nn-i+l  K(cB,n^  *  1  /  1  n-i+1  R^CB  +  • 

From  (A3)  we  can  easily  obtain 

J0  P,VCB>  *  11  jQ  P{ZB<CB  +  »  *  11 
*  P|W  ■  J>  -  "p(ZB(c8-H)  ■  JJ 

which  implies,  when  taking  reciprocals, 

P{ZB(CB)  *  j  1  ZB(CB)  ~i}  ^  P{ZB(CB  +  X>  “  J  I  ZB(cB  +  5  J>  • 

which  in  turn  implies,  by  subtracting  both  sides  from  one. 
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P^ZB^CB)  =  I  2B(cB)  "  N  piZB(cB+l>  v  J-l  I  .1*  . 

1  <  j  <  Hty  . 


P{W  <  J-l)  =  P{ZB(cB)  5  j-1  |  ZB(cfi)  <  M+y} 

M+y 

’  j”  P{W  s  I  W  i  o 

M+y 

-  ifj  P{ZB(CB+1)  -  i-l  I  ZB(V1)  £  i} 

-  p{zB (cB+i)  <  j-l}  . 

Hence , 

P{2B(cB)  -  j}  i  p(ZB(cB+1)  >  j}  . 


Ordering  the  ratios  as  in  (A2)  actually  is  a  sufficient  condition 
in  general  for  stochastic  ordering  [see  WHITT  (1980)]. 
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To  cope  with  the  expanding  technology,  our  society  must 
be  assured  of  a  continuing  supply  of  rigorously  trained 
and  educated  engineers.  The  School  of  Engineering  and 
Applied  Science  is  completely  committed  to  this  ob¬ 
jective. 
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